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ABSTRACT 

We construct exact solutions to Einstein equations which represent relativistic disks 
imbedded in an expanding FRW Universe. It is shown that the expansion influences the 
kinematical characteristics of the disks such as rotational curves, surface mass density, etc. The 
effects of the expansion are exemplified with non-static generalizations of Kuzmin-Curzon and 
Schwarzschild disks. 


Subject headings: celestial mechanics, stellar dynamics — cosmology: large scale structure of the 
Universe — galaxies: kinematics and dynamics — ISM: kinematics and dynamics — relativity 


1. Introduction 


Recently, there has been a renewed interest in the study of relativistic disks (Bicak, Lynden-Bell & 
Katz, 1993; Bicak & Ledvinka, 1993; & Bicak, Lynden-Bell & Pichon, 1993), stemming from a desire to 
address the possible observational evidence of giant black holes sustained by the surrounding disks. The 
gravitational fields of galactic disks may be accurately modeled by the Newtonian potential theory provided 
the thickness of the disk is negligible compared with the typical size of their halos. In more sound situations, 
however, when gravity is strong enough one must turn to General Relativity. A typical situation, where it 
is believed that Einstein theory may contribute is met when one considers the accretion disks around the 
central black holes in quasars. Disks may be also used to model sheet-like structures, so their study may 
shed some light on the understanding of the large scale inhomogeneities present in the Universe (Lemos & 
Ventura 1994). 


In General Relativity, solutions representing thin disks may be constructed starting with static axially 
symmetric vacuum Weyl metrics. A discontinuity is then introduced in the first derivative of the metric 
across the azimuthal plane z = 0, which in turn induces distribution-like terms in the Ricci tensor. 
Everywhere outside the z = 0 plane the solution is vacuum with a non-vanishing surface mass density 
concentrated on this plane. Since Weyl solutions are static, they should represent static disks. To allow for 
rotation in General Relativity one must drop the assumption of the orthogonality of the two Killing vector 
fields associated with Weyl geometry, thus complicating the problem by introducing the dragging of inertial 
frames. 


Fortunately, Morgan and Morgan ( |l969a| ; [1969b ; 1970 ) had the very smart idea of interpreting the 
“static disks” as made out of two equal streams of collisionless particles circulating in opposite directions. 
Consequently, the total angular momentum vanishes and the system may be described by a Weyl line 
element. Interestingly enough, there exists observational evidence of some galaxies with two counterrotating 
stellar components around their center (Rix, Franx, Fisher & Illingworth 1992; Merrifield & Kuijken 1994; 
Bertola et al. 1996; Kuijken, Fisher & Merrifield 1990). 
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Another interpretational problem with the solutions described by the Weyl metric arises due to the 
fact that the disks are infinite in their extent. One may argue, however, that infinite disks model the inner 
portions of galaxies or accretion disks. Once the interpretational problems are dealt with, it is not difficult 
to model relativistic disks and study their physical properties such as velocity profiles, surface mass density, 
redshifts and so on. 


Evans and de Zeeuw ( 1992| ) showed that it is possible to analyze any classical axially symmetric disk 
into a linear superposition of so-called Kuzmin disks ( Kuzmin 1956| ). Using this result, Bicak, Lynden-Bell 
and Katz ( 1993| , hereafter BLK) constructed families of counterrotating disk space-times. Their work 
was further extended by Bicak, Lynden-Bell and Pichon (1993, hereafter BLP) and more recently by 
Lynden-Bell and Pichon (1996). On the other hand, Lemos and Letelier (1993; 1994; 1996) using the well 
known techniques of superposing different Weyl solutions, have constructed space-times describing disks 
surrounding Schwarzschild black holes. 


The main purpose of this paper is to consider, in the framework of exact solutions in General Relativity, 
the effect of the cosmological expansion on the disk dynamics on one hand, and the influence of disk-like 
large-scale structures on the cosmological model’s kinematics on the other. It is known (Noerdlinger & 
Petrosian 1971 ) that the effects of the expansion on bounded gravitational systems is proportional to 
the ratio p/p s , where p is the background energy density and p s is the mean rest-density of the system. 
Therefore, for galactic and accretion disks fueling the quasars these effects are not expected to be too 
significant, though the overall expansion may change the disk dynamics. For disk-like large-scale structures, 
however, the expansion may turn out to be of particular importance. 


In §2 we propose an algorithm to generate solutions to Einstein equations which may be interpreted 
as inhomogeneities imbedded in a spatially flat isotropic Universe. The technique uses a scalar field 
generalization in a cosmological setting of static Weyl solutions. The scalar field can be split into two parts, 
one homogenous, which may be hydrodynamically interpreted as the velocity potential for an adiabatic 
perfect fluid acting as the source of a FRW expansion; and another one highly irregular which may be 
thought of as a local disk inhomogeneity. In §3 we apply the generating procedure of the previous section 
to the general relativistic Kuzmin disks. We then analize the energy-momentum tensor for the obtained 
solutions and interpret them as disks imbedded in a FRW background. §4 is devoted to study the dynamics 
of the disks with the emphasis on the effects of the expansion on pertinent dynamical and kinematical 
quantities. We conclude the paper with the discussion and future prospects. 


2. Exact non-static solutions with axial symmetry 
2.1. The generating algorithm 


In this Section we obtain non-static solutions to Einstein field equations representing compact objects 
in a cosmological setting. We start with Weyl’s line element which can be written as: 


ds 2 = -e 2 ^ p ’ z) dt 2 + e~ 2uM p 2 dcl) 2 + e~ 2uM+2<:{p ’ z) {dp 2 + dz 2 ). 


(1) 


Throughout the paper we use G = c = 1. It is well-known that vacuum static axially symmetric fields in 


General Relativity can be generated starting with a Newtonian potential ( Weyl 1917 ; Levi-Civita 1919a ; 
Levi-Civita 1919b). For the space-times obtained in this way the metric function v in equation (Q) may be 
taken to be any classical solution of the Laplace equation in cylindrical coordinates, and the other metric 
function £ is obtained by a quadrature. 
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New non-static metrics may be generated in two stages. First, starting from a vacuum solution of the 
Einstein equations for a metric given by equation (if]), we construct a new static solution with a minimally 
coupled massless scalar field ip as a source. For these space-times the energy-momentum tensor takes the 
form: 


Tab = Ip,alp,b - 9ab Ip,c4> , 

(2) 

and the set of Einstein equations is 

b*ZZ + Vpp T P Up — 0 , 

(3) 

Czz + C pp ~ P Cp T 2 Up = —ipp , 

(4) 

Czz + Cpp + P 1 Cp + 2 U 2 = — Ipz 2 , 

(5) 

P -1 Cz -2 u p u z = -Ipzipp , 

(6) 

along with the Klein-Gordon equation 

Ipzz +1ppp + P 1 lpp = 0 . 

(7) 


It is easy to see that if u a and ( 0 solve the vacuum Einstein equations, a solution to equations (||-^) is then 
given by 

v = u Q + C log p , (8) 

C = B Co + E u a + F log p , (9) 

ip = Au a + D log p , (10) 

where the constants are subject to the following constraints: 

2C + AD = E, (11) 

A 2 + 2 = 2 B, (12) 

D 2 + 2 C 2 = 2 F. (13) 

For the purposes of this work we will discard the logarithmic terms by setting the constants C, F and D to 
zero, ensuring thus asymptotic flatness at spatial infinity in a static case. 

Once the solution with a massless scalar field metric has been constructed, we transform it into a 
non-static solution with a self-interacting scalar field with an exponential potential. This is accomplished 
by using an algorithm originally due to Fonarev ( 1995| ); see as well the generalization by Feinstein, Ibanez 
& Lazkoz (1995|). The new solution in a synchronous system of coordinates reads 


ds 2 = -e 2uM dt 2 + R 2 {t)e~ 2 ^ p ' z) (p 2 d(P 2 + e 2 ^ p ' z \dp 2 + dz 2 )) , 

where the “scale factor” is R(t) = f 2 / fe2 and the new metric functions are 

u = z'o and £ = (1 + 2/A: 2 ) Co , 


(14) 


(15) 


the constant k being the slope of the potential V = Ae and A = (12 — 2 k 2 )/k 4 . The line element in 
equation ( 0 ) is a solution of the Einstein field equations with the energy momentum tensor given by 


Tab = i>,a1p,b ^ gab QV’.iV’’* + k . 


(16) 
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Furthermore, the new scalar field rl>(t , p , z ) splits into a homogeneous and an inhomogeneous part, 
ip(t, p, z) = 4>h + 4>inh with ip h and ip inh given by 

2 2 

i’h = -j; logt and ipinh = - Vo- (17) 

Note that in the particular case k 2 = 6 the potential term vanishes, and one is left with a massless scalar 
field. 


2.2. Interpretation of the new solutions 


As long as the seed static metric is asymptotically flat, the newly generated non-static solution is 
asymptotically homogeneous and isotropic, and at spatial infinity the geometry is represented by the FRW 
metric 

ds 2 = — dt 2 + ( p 2 dq i > 2 + dp 2 + dz 2 ) . (18) 

We now identify the homogeneous part of the scalar field with the velocity potential of an irrotational 
perfect fluid which is the source of the metric in equation ( |l 8 |) . The energy-momentum tensor is then 

(19) 


Tab RW = (P + k)Ua'Ub +pgab, 
where the four-velocity of the fluid is given by 

_ i>h,aS? 

V(—0M h,t ) ’ 

so that the pressure p and the energy density g become 

M = ~^h,tfp h ’ t + and p = tl> h>t i’ h ^ - At~ 2 . 


Substituting the expression for the scalar field we readily get 

4 k 2 - 12 


P = 


and g = 


12 


t 2 fc 4 r t 2 k A 

It is straightforward to see that equation (|2^) defines the barotropic equation of state 

k 2 — 3 


p = 7 g and 7 = 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


Note, however, that the scalar field has also an inhomogeneous component. Although the 
inhomogeneities steadily dilute as the spatial distance increases, they cannot be regarded to be negligible in 
the intermediate regions. 

It is well-known that as long as the gradient of the scalar field is timelike the energy momentum tensor 
can be interpreted in terms of a perfect fluid. The asymptotic homogeneity of the solution ensures that 
at any time the gradient of the scalar field is timelike outside a certain closed spatial region, depending 
basically on the mass and compactness of the disk, as well as on the strength of the scalar field. Outside 
this region the perfect fluid interpretation holds, and at spatial infinity the fluid becomes homogeneous and 
isotropic as mentioned above. 
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On the other hand, the evolution of the local inhomogeneities with time may also be deduced by 
studying the gradient of the scalar field: 


^,oV'’ C 


Iph, c4’ h ’ C + ifrinh, c1p mh ’ C , 


(24) 


where 


and 


^h,ci> h ' c = , t 2 e“ 
3y + 3 


ink. cP 


37 + 3 


t + v 2 z ). 


(25) 

(26) 


We have checked numerically that for the models with accelerated expansion 7 < —1/3 the outer 
region, where the homogeneous term of the scalar field dominates the one representing the spatial gradients, 
grows with time. On the other hand, the highly irregular inner region shrinks to zero. This behaviour is 
similar to that one found in the inflationary scenario, in the sense that the small-scale inhomogenities are 
washed away due to the accelerated expansion. 

The highly inhomogenous region surrounding the disk plane where the perfect fluid interpretation is 
not valid may be thought of as playing the role of a cushion between the disk and the cosmological fluid, 
representing a transition region between both regimes. 

We suggest, therefore, to interpret the solutions of the type given by equation (0 as describing an 
inhomogeneity embedded in an expanding Universe on the grounds that on large scales one recovers the 
Friedmannian behaviour while the local structure for small distances is governed by the seed metric. This 
interpretation will be re-enforced in the light of the results of following sections. 


3. Relativistic expanding thin disks 

We now apply the described algorithm to the well-known family of axially symmetric solutions 
representing infinitesimally thin disks in a relativistic context and construct their dynamical counterparts. 
We are mostly interested in the influence of the expansion on the kinematical quantities that characterize 
the disks. It is expected that for an expanding Universe the mass-energy density on the disk plane will 
steadily decrease as time goes by. However, the effect on other quantities such as velocity or angular 
momentum is not intuitively foreseeable. In addition, we believe that the conclusions reached in this work 
can bring a new perspective onto the problem of embedding irregular sources into standard cosmologies. 


3.1. Generation of thin disk families 


The disk configurations we are dealing with were firstly studied by Morgan & Morgan ( 1969a , 1969b , 
1970|) and describe counterrotating disks with the same number of leftwise and rightwise rotating particles. 


As was pointed out in the previous section one can start with the Newtonian potential produced by a thin 
disk to obtain its relativistic version. 


For the sake of simplicity, unless otherwise stated, we will focus the analysis on two families of disks, 
the Kuzmin-Curzon disks and the generalized Schwarzchilcl disks (BLK, 1993). The Kuzmin-Curzon disks 
may be used as building blocks for infinite families of more complicated disk-like solutions; one should just 
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superpose elementary Kuzmin-Curzon disks with different compacticities weighted by a function W(b), as 
shown by Evans & de Zeeuw (1992). Their work on Newtonian axisymmetric potential-density pairs was 
extended to a relativistic case by BLK, who gave the general form of the metric for a superposition of 
Kuzmin-Curzon disks. 


BLP gave a step further in the study of relativistic disks by constructing the more involved families 


of relativistic versions of the Kuzmin-Toomre (Toomre 1963; Nagai fe Miyamoto 1976; Evans & de Zeeuw 


1992) and Kalnajs-Mestel disks ( Mestel 1963 ; Kalnajs 1976 ). The latter family of solutions is particularly 
interesting because they have long and flat rotation curves for certain values of the parameters, and include 
the so-called generalized Schwarzschild disks as their lowest order representatives. 

The potential corresponding to a Kuzmin-Curzon disk is obtained by introducing a discontinuity in 
the original Curzon metric ( |Curzon 1924 ; Chazy 1924 ) using the transformation z —> \z\ + b. This way of 
introducing the discontinuity is equivalent to placing two mirror particles of mass M at a distance b below 
and above the 2 = 0 plane on the 2 axis. The Newtonian potential for the Kuzmin-Curzon disk is 


K-C 


M 


[ P 2 + (\z\+bn- 


(27) 


and following Weyl (1917), one identifies it with the metric function v a in equation This metric function 
remains a solution of the Laplace equation outside the 2 = 0 plane, but now due to the discontinuity, there 
appears to be a non-vanishing surface mass density. Integration of the remaining metric function yields 


fK-C _ 


M 2 p 2 


[p 2 + (M + &)T 

Correspondingly, as shown by BLK, the classical surface mass density £(p) = 2v 0 | 


^K—C 


4Mb 
(p 2 + b 2 ) z 


(28) 


(29) 


Here M is the total mass of the disk as measured from the infinity, and b is a parameter measuring the 
compactness of the disk. 

The analogue of a generalized Schwarzschild disk is constructed in the classical Kuzmin’s picture by 
substituting the two mirror particles by two rods with constant line density. Then, the general relativistic 
solution is given by 

Pmin + | z\ T b rn/ i 


v?~ S = 


M 

bmn.x b 11 


■log 


and 


pG-S 
So 


2 M 2 


■log 


(Pn 


Pmax + |~| + bmax 
T Pmax) (bmax ^mm) 


(b - b ■ ) 2 

\ u max u min) 

In this case, the classical surface mass density D(p) is 


4 Pmin Pn 


(30) 

(31) 


e g-s 


4 M 


(bmax bmin) (p 2 T b 2 )‘ 


(32) 


where p^ in = p 2 + (\z\ + b min ) 2 and p^ irl = p 2 + {\z\ + b min ) 2 , b max - b min being the parameter that 
measures the compactness of the disk. 
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3.2. Disks in an expanding FRW Universe 

We start by looking somewhat more carefully at the stress-energy tensor for the solutions given by 
equation © with the metric functions corresponding to a generic thin counterrotating disk. 

In order to deal with the discontinuities in the metric across the z = 0 plane, we consider the metric 
functions in the sense of distributions and introduce a new variable £ = \z\. Generically, the metric 
coefficients will have square-integrable weak derivatives, and so the usual formula for the Ricci tensor can 
be interpreted in the sense of distributions. The stress-energy tensor reads (cf. as well Chamorro, Gregory 
& Stewart, 1987): 


rj ip rpZ 

p ^ Z 

= 2R(ty 2 (p-\, p + ^-u 2 p )e- 2 ^, 

(33) 

'TP r p z 

- 1 - p ' ± z 

= —2 R(t)~ 2 [r 2 + 2RRj e~ 2 v , 

(34) 

T$+T} 

= — 2 R(t) 2 [y,p P + + 2 8(z) + p 1 i' tP — 

- 2 S(z) C, £ - v 2 p - 4+) e" 2 c+2 u 



- 2 R(t )~ 2 | 

[2 R 2 + RR^J e ~ 2l/ , 

(35) 

T t _ T * 

— —2R(t) 2 [y,p P + + 2 8(z) v t £ + p 1 R,p) e 2 ?+ 2!y 



- 2 R(t)~ 2 | 

[r 2 -RR) e~ 2v , 

(36) 

- z 

= —2 R(t) 2 [8(z)~ 0(-2)] (2i/ iP ^ - p 1 C,c) e 2C+2y , 

(37) 

T p 

= 2 R{t)~^Rv p e -2< U 2ly , 

(38) 

rpZ 

1 t 

= 2 R(t)~ 3 [8(z) - 8(—z)} e~ 2<:+2 r 

(39) 


The energy-momentum tensor given by equations ( [l3f]39| ) splits into a regular and a singular part 
proportional to the functional 5(z), which is interpreted as a thin counterrotating disk. We still separate 
the regular part into two bits: a highly inhomogeneous part vanishing at infinity that we interpret as an 
interaction term between the matter composing the disk and the cosmic fluid, and another part representing 
at infinity an isotropic perfect fluid. Phenomenologically we write 

rj-ib _ rj-ibdisk _|_ rj-ibint _|_ 2 i"j-ibFRW ( 40 ) 

where 

T} FRW = 3 R{t)~ 2 R 2 , (41) 

and 

T gFRW = T pFRW = t zFRW = _ 2 R ( t y 2 ^2 + 2 R ^j _ (42) 

The identification of the other terms is straightforward. The term we have referred to as X), ,sfc is precisely 
the one which is discontinuous across the z = 0 plane and which therefore gives rise to the non vanishing 
mass density. One has to be careful when interpreting equation (^), for it has a purely phenomenological 
character. Obviously none of the three terms in the r.h.s. of equation ( [f0| ) can be regarded as a true 
energy-momentum tensor, since they do not satisfy the energy conservation equation T ab = 0 separately. 
Note, that this decomposition into three terms is similar to that proposed in §2.2 based on the character of 
the gradient of the scalar field. Both approaches qualitatively lead to the same physical description. 


4. Dynamics of expanding disks 


In this Section we study the kinematical quantities of relativistic disks and discuss the influence of 
the expansion. We explicitly prove that the mass-energy density decreases with time, and also address the 
question of the effect of the expansion on the flattening of the rotation curves. Special attention will be paid 
to the comparison between a static disk and its counterpart in a pressure-free Universe. Our motivation is 
that the Universe, as observed at present, has negligible pressure. 

We now look at the surface mass density and streaming velocities of the particles on the disk. 

The r b surface components are obtained by integration across the disk of the T' b components of the 
energy-momentum tensor, 

r b a = I T b R(t)e^dz. (43) 


As shown by BLK, the surface rest mass density in a fixed reference system reads ao = 2cr p (l — v 2 ) -1 / 2 , 
while the surface mass density in the same reference system is <r = 2 cr p (l — u 2 ) -1 = <ro(l — u 2 ) -1 / 2 , where v 
is the rotational velocity and a p is the proper rest mass density of one stream. 

We first calculate the surface mass density cr(p,t ) on the disk plane by integrating the T\ term across 
the disk. Using equations (|3|,|3(]) one gets 

/* 0 + 

o"(p, t) = —t\ = - Tt i?(f)e (c_1/) dz = i?(f) _ 1 e (l/_c) (4 - 2 C,£)|«=o • (44) 

Jo- 

Substituting (y as calculated from equation (£J) we obtain after some algebra 


<j{p,t) =4R(t) C) (l-p/(7 )Up)I«=o, 


(45) 


where we denote /(y) = (3 7 + 3)/(3 7 + 5) and 7 is the adiabatic index of the perfect fluid. Note that the 
static case will be recovered in the limit 7 —> 00 , /(y) —* 1. The particular case of a dust filled Universe 
corresponds to 7 = 0 and /( 7 ) =3/5. 

By integrating the term obtained from equations (^.@ one gets 


V (T = Tl = 


0 + 


0 - 


T/R(t)e^dz = 2 R(t)~ 1 e^f(' Y ) Q |*=o , 


(46) 


or alternatively 


v 2 a = 4 R(t) V" C) p/(y)Up^lf=o • (47) 

Then, from equations ( |45| , [47| ) one readily obtains the expression for the square of the streaming velocities 

P/(7) Up 


v 2 = 


1 ~ Pfil) U/ 


(48) 


£=o 


In order to emphasize the difference in the flattening of the rotation curves in the static and expanding 
cases we also evaluate the derivative of the velocity in the radial direction obtaining 


/(t) 1/2 ^. 




2[Un(i-p/(y) U p )] 1/2 


(49) 


«=o 
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Since the classical mass density £(r) = 4 i/^|j = o of the disks is positive everywhere, in order to ensure 
the positivity of the surface energy density a we must impose p/(y) ^. p \^=o < 1, as deduced from equation 
([f5|). This is equivalent to the fulfillment of the weak energy condition (cf. BLK). For the particular case of 
an expanding Kuzmin-Curzon disk this restriction reduces to 


Mf( 7 ) \/27 

b 2 

whereas for a generalized Schwarzschild disk it becomes 


M /(y) 


(Pmax bmin ) 


< 


[1 (bmin/bmax)~\ 2 
[1 {bmin/bmax) 3 ] 2 


(50) 


(51) 


therefore, the expanding disks may be denser than their static counterparts. 

Moreover, if p/(y) v tP |^ = o <1/2 the dominant energy-condition holds and the rotation velocities are 
subluminal. It is clear that as long as the dominant energy condition holds the positivity of the mass 
density is also ensured. 

To compare the streaming velocities of a static disk and any of its time-evolving counterparts, one may 
look at the expression of their quotient, which reads 


v 2 (-y = oo) 


/( 7 ) 


1 ~ pv, P 
1 — P f{l) v ,p 


«=o 


(52) 


Since /(y) < 1, we conclude that for expanding disks the velocities are typically everywhere lower than 
those for a static one (see Figure 1). 

To compare our results with those obtained by BLK we represent the kinematical quantities as 
functions of the circumferential radius p c = pe - ", where 2irp c is the physical circumference of a circle with 
radius p = constant. Note that for the solutions studied in this paper the comoving radius does not depend 
on the adiabatic index y. 


Yet, more interesting conclusions are obtained from the quotient between the velocity derivatives in 
the radial direction, 


V 'Pb) = f l/2, ) ( l ~PV,P X' 2 

V A y = oc) 7) \i - pfA)v, P J 


(53) 


From equation (|5^) it can be learnt that the slope of the rotation curves for expanding disks is less 
pronounced than for static ones. It can be concluded that the expansion accentuates the flattening of the 
rotation curves (see Figure 2). 


Another question to dwell on is the effect of the expansion on the mass-energy density on the disk 
plane. Similarly, as done above with the velocities and their radial derivatives, we look at the quotient of 
the mass energy density of an evolving and a static disk: 


= ( ie C0)-3W3 

<r(y = 00 ) 1 ~ pv, P 


(54) 


Since Co is negative everywhere on the disk, we conclude that expanding disks are denser than static ones, 
with additional mass concentrated outside the center, whereas the central density is not affected by the 
expansion. Time-evolution of this quantity crucially depends on the adiabatic index y. To be more specific, 
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the mass-density dilution rate increases, keeps constant or decreases with time, if the expansion parameter 
7 is smaller, equal or larger than —1/3 respectively. 

It is also interesting to examine the surface rest mass density of the rotating particles as measured in a 
frame attached to fixed axes. This quantity is given by 

cr 0 = a (l - ii 2 ) 1/2 . (55) 

Substituting the expression for a and v in equations (^5j,^8)) we get 

cr 0 (p, t)= 4 e (l/_c) [1 - p f{j) v p \ i [1 - 2 p /( 7 ) u p ] 3 | e=0 . (56) 


The latter quantity is everywhere positive provided the dominant energy condition holds. Taking the 
previous results into account it is straightforward to see that the influence of the expansion on the surface 
rest mass density a o is the same as for the mass density a. 

We have further defined a new function Act/ct = [<7(7) — <7(7 = oo)]/ct(7) which gives the contrast 
in surface mass density between an expanding and a static disk. Figures 3a and 3b represent the 
contrast function Act/ct and Figures 4a and 4b give their rest counterparts Acto/cto for both a generalized 
Schwarzschild disk and a Kuzmin-Curzon disk. It can be learnt from the figures that the expanding disks 
with the same mass parameter M are much denser than their static counterparts. 


Finally, it is worth looking at the specific angular momentum of the counterrotating components of the 
disk. Bearing in mind that the specific angular momentum of a particle with rest mass m rotating at radius 
r is defined as 

j = (P4><p/‘m) = g^^dfi/dX, (57) 


where A is proper time, then 


or alternatively 


j(p,t) = R(t ) 


pv e 


(1 — v 2 )^ ^ 


j(p,t) = R(t ) 


(P 3 f( i)v,p)'e~ v 


(1-2^/(7) Fp) 2 


(58) 


(59) 


e=o 


Similarly to what happened with the rotation velocities, a smaller 7 gives a smaller angular momentum, 
which increases with time. This behaviour is presented in Figures 5a and 5b where we compare the 
specific angular momentum in an expanding and a static case for a generalized Schwarzschild disk and a 
Kuzmin-Curzon disk. 


5. Conclusions and outlook 

We have constructed exact solutions to Einstein equations which we interpret as representing relativistic 
disks in a cosmological setting. A self interacting scalar field serves as a source of the global expansion, yet 
locally it defines a disk-like structure across the z = 0 plane. Far away from this plane the gradient of the 
scalar field may be identified with the velocity potential of an irrotational perfect fluid with an adiabatic 
equation of state p = 7 p. Near the azimuthal plane the spatial gradients of the scalar field dominate over 
the kinetic part and the geometry is highly inhomogeneous. No sharp-cut transition region exists in between 
these two different regimes, however, this probably corresponds to a physically more realistic situation than 
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a surface matching between different solutions of gravitational field equations. Moreover, to the best of our 
knowledge, there do not exist such solutions in the case of axial symmetry and an external FRW geometry. 


Some particular cases of the solutions we have obtained ( k 2 = 6) may be re-interpreted as disks in 
the Brans-Dicke theory. Indeed, in this case the scalar field may be considered as massless (the potential 
term vanishes) and a simple conformal transformation ds 2 —* \^“ +3 / 2 ds 2 transforms the solution into a 

Brans-Dicke frame (Tabensky and Taub 1973). The physical interpretation of the solutions then is quite 
different: these represent relativistic disks in a theory where the gravitational constant varies in space and 
time. We will not dwell more about this point here, but just mention that the analysis of our paper applies 
equally well to Brans-Dicke disks. 


Once the solutions are interpreted as local inhomogeneities in a model Universe the way is cleared 
to see the effects produced by expansion. We have found that depending on the rate of expansion the 
inhomogeneities occupy larger or smaller regions of the Universe. More specifically, for the accelerated 
expansion we find that inhomogeneities disappear with time, while in a deccelerated model their growth is 
unbounded. 


The effect of the expansion on the kinematical characteristics of the disks was studied as well. We have 
shown that expansion changes in principle the fall-off of the rotational curves, the angular momentum and 
the surface mass density of the disks. We have also compared the characteristics of the static disks with 
those in a dust filled Universe. 

Although we have concentrated our study on disk-like objects, it is remarkable that the generating 
technique used to obtain solutions is also adequate to study other type of sources of astrophysical interest 
such as cosmic strings, walls, spherical shells, etc. 
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Fig. 1.— Streaming velocity v as function of the circumferential radius p c for a generalized Schwarzschild 
disk with M = 100, b m i n = 1 and b max = 200 (a); and for a Kuzmin-Curzon disk with M = 1 and b = 4 (b). 
The continuous lines correspond to expanding disks in a dust filled FRW Universe,whereas the dashed lines 
correspond to their static counterparts. 
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Fig. 2a Fig. 2b 


Fig. 2.— Radial derivative of the streaming velocity v p as function of the circumferential radius p c for a 
generalized Schwarzschild disk with M = 100, b m i n = 1 and b max = 200 (a); and a Kuzmin-Curzon disk with 
M = 1 and 6 = 4 (b).The continuous lines corresponds to expanding disks in a dust filled FRW Universe, 
whereas the dashed lines correspond to their static counterparts. 






Fig. 3a 


Fig. 3b 


Fig. 3.— Surface energy density contrast Acr/er at t = 1 as function of the circumferential radius p c for a 
generalized Schwarzschild disk with M = 100, b m i n = 1 and b max = 200 (a); and for Kuzmin-Curzon disk 
with M = 1 and b = 4 (b). 
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Fig. 4b 


Fig. 4.— Surface rest energy density contrast A cto/cto at t = 1 as function of the circumferential radius p c 
for a generalized Schwarzschild disk with M = 100, b m i n = 1 and b max = 200 (a); and a Kuzmin-Curzon 
disk with M = 1 and 6 = 4 (b). 
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Fig. 5a Fig. 5b 

Fig. 5.— Specific angular momentum j at t = 1 as function of the circumferential radius p c for a generalized 
Schwarzschild disk with M = 100, b m i n = 1 and b max = 200 (a); and for a Kuzmin-Curzon disk with M = 1 
and b = 4 (b). The continuous lines correspond to disks in a dust filled FRW Universe, whereas the dashed 
lines correspond to their static counterparts. 





